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ABSTRACT 


This  MTR  contains  a  summary  of  some  of  the  more  important  theoretical 
results  related  to  the  performance  of  a  digital  matched  filter  receiver  de¬ 
signed  to  receive  a  biphase  modulated  signal  in  the  presence  of  interference. 
Both  the  synchronization  and  data  performance  of  the  receiver  are  evaluated 
when  the  signal  is  received  in  the  presence  of  a)  Gaussian  noise,  b)  Random 
phase  sinusoidal  interference,  c)  Fixed  phase  sinusoidal  interference.  Specific 
attention  is  paid  to  the  effect  of  processing  the  received  signal  using  one  bit 
quantization. 
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SECTION  I 


INTRODUCTION 

This  MTR  is  a  summary  of  some  of  the  more  important  results  re¬ 
lated  to  the  performance  of  the  PN  modem  being  built  by  D-8  of  the 
MITRE  Corporation.  MTR-2446  by  R.  Haggarty  will  contain  the  detailed 
circuit  design  of  the  modem  together  with  details  of  the  RF  transmitter 
and  receiver  units  and  laboratory  test  data  on  the  modem  performance. 

In  order  to  increase  the  utility  of  the  work  in  this  MTR,  where 
possible,  the  results  obtained  have  been  stated  in  terms  of  the  important 
parameters  of  a  more  general  model  of  the  modem  than  that  actually  being 
built.  Since  the  modem  was  designed  for  use  in  a  tactical  environment 
emphasis  is  placed  on  the  anti-jam  performance. 

In  this  communication  system  information  in  binary  form  is  trans¬ 
mitted  by  sending  one  of  two  uncorrelated  signals.  The  signals  are 
derived  by  choosing  one  of  two  pseudo-random  sequences,  which  is  used 
as  a  rectangular  waveform  to  directly  phase  modulate  a  sine  wave  carrier. 
The  optimum  receiver  consists  of  two  matched  filters,  each  matched  to 
one  of  the  two  possible  transmitted  signals,  followed  by  a  level  com¬ 
parison  circuit.  This  implementation  assumes  that  the  receiver  has  de¬ 
rived  the  received  phase  of  the  RF  carrier.  If  it  does  not,  the  matched 
filter  for  each  of  the  two  possible  signals  is  replaced  by  two  filters 
each  matched  to  the  same  signal  but  90°  out  of  carrier  phase  with  one 
another.  The  outputs  of  each  pair  of  quadrature  matched  filters  are 
squared  and  added  before  entering  the  level  comparison  circuit. 

A  matched  filter  for  the  above  mentioned  signal  is  commonly  im¬ 
plemented  by  a  mixer  to  convert  the  RF  signal  to  baseband  followed  by 
an  analog  tapped  delay  line,  with  plus-and-minus-one  tap  gains  corres¬ 
ponding  to  the  appropriate  pseudo-random  sequence,  and  connected  to  a 
single  adder.  The  tapped  delay  line  can  be  replaced  by  an  analog-to- 
digital  converter  followed  by  a  Q-level  quantizer  and  a  shift  register 
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with  taps  again  matched  to  the  appropriate  pseudo-random  sequence.  If 
the  number  of  levels  Q  used  by  the  quantizer  is  large  the  performance 
of  the  receiver  approaches  that  of  an  analog  delay  line  implementation. 
Both  the  analog  delay  line  and  the  required  amount  of  storage  for  a 
digital  implementation  using  a  large  number  of  voltage  levels  Q,  tend 
to  be  expensive.  The  simplest  and  cheapest  implementation  is  do  to  it 
digitally  using  only  two  quantization  levels  (i.e.,  Q  =  2),  however, 
one  expects  to  pay  a  price  in  performance  for  this  simplification. 

The  majority  of  the  results  given  in  this  MTR  are  concerned  with 
the  performance  both  in  the  ^nchronization  and  data  modes  when  the 
received  signal  is  processed  using  one  bit  quantization. 
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SECTION  II 


DESCRIPTION  OF  THE  MATCHED  FILTER  RECEIVER  MODEL 

2 . 1  Basic  Matched  Filter 

Figure  1  shows  a  block  diagram  of  the  basic  digital  matched 
filter.  The  received  bi-phase  modulated  sequence  is  separated  into 
inphase  and  quadrature  baseband  components  by  the  two  mixers  shown  in 
the  diagram.  Both  the  inphase  and  quadrature  components  are  low  pass 
filtered  before  being  hard  limited  by  the  1  bit  A/D  converters.  The 
output  of  the  two  one  bit  A/D  converters  are  sampled  at  the  PN  rate 
and  the  "ones"  and  "zeros"  (corresponding  to  plus  and  minus  voltages) 
are  shifted  into  the  appropriate  shift  register.  A  copy  of  the  trans¬ 
mitted  PN  sequence  is  stored  in  comparison  shift  registers  and  the 
contents  of  the  I  and  Q  channel  shift  registers  are  continually  com¬ 
pared  against  the  contents  of  the  comparison  registers.  The  sum  of 
agreements  minus  disagreements  in  both  the  I  and  Q  channels  (i.e.,  V ^ 
and  V^)  are  continuously  calculated,  squared  and  added  and  the  square 
root  of  the  sum  Z  calculated.  When  the  received  sequence  is  correctly 
aligned  in  the  I  and  Q  shift  registers  there  is  a  large  output  from 
the  square  root  circuit,  i.e.,  the  received  sequence  is  highly  correlated 
with  the  stored  replica. 

2 . 2  Binary  Receiver  Using  Two  Matched  Filters 

For  the  transmission  of  binary  information  two  orthogonal  PN 
sequences  would  be  selected,  one  corresponding  to  a  binary  "one",  the 
other  to  a  "zero".  The  receiver  would  contain  two  digital  matched 
filters  as  shown  in  Figure  2,  one  filter  being  matched  to  the  "zero" 
sequence  and  the  other  matched  to  the  "one"  sequence.  Assuming  that 
the  receiver  is  synchronized  the  outputs  Z  and  Z^  from  the  two  filters 
would  be  compared  at  the  end  of  a  received  information  bit  and  a  decision 
as  to  the  reception  of  a  "one"  or  "zero"  would  be  based  on  the  larger 
of  the  two  outputs  Z  or  Z^. 
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Figure  I  BASIC  DIGITAL  MATCHED  FILTER 
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Figure  2  BINARY  RECEIVER  USING  2  DIGITAL  MATCHED  FILTERS 


SECTION  III 


PROBABILITY  DENSITY  FUNCTION  FOR  THE  OUTPUT  VOLTAGE  OF  THE  DMF 


In  order  to  compute  such  performance  parameters  as  the  bit  error 
probability,  the  probability  of  falsely  synchronizing  and  the  probability 
of  missing  sync  we  first  need  the  probability  density  function  for  the 
output  voltage  from  the  digital  matched  filter  (DMF).  In  this  section, 
we  give  this  probability  density  function  when  the  bi-phase  modulated 
signal  is  received  together  with  additive 

(a)  white  Gaussian  noise 

(b)  random  phase  sinusoidal  interference 

(c)  fixed  phase  sinusoidal  interference 

The  actual  derivation  of  the  probability  density  functions  are  given 
in  Appendices  A,  B  and  C. 

In  order  to  simplify  the  analysis  of  the  DMF  performance  the 
following  simplifying  assumptions  were  made  in  the  appendices.  They 
are, 

(1)  that  the  low  pass  filter  in  each  channel  is  wide  enough 

to  allow  the  pulse  amplitude  to  build  up  to  full  amplitude 

* 

within  a  chip  time.  In  the  actual  design  of  the  filter  there 
will  be  the  usual  conflict  of  requiring  the  filter  to  be  wide 
to  allow  fairly  rapid  build  up  to  full  pulse  amplitude  but  at 
the  same  time  keeping  it  as  narrow  as  possible  to  reject  as 
much  noise  as  possible. 

(2)  that  sampling  of  the  I  and  Q  channel  A/D  converters  takes 
place  towards  the  end  of  each  chip,  that  is  when  the  low  pass 
filter  outputs  have  reached  maximum  amplitude.  To  ensure  that 


A  chip  time  is  the  duration  of  the  basic  signaling  element  and  consists 
of  the  transmitted  carrier  in  one  of  two  possible  phases  (0°  or  180°). 
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this  is  achieved  it  is  usual  to  duplicate  the  matched  filter 
assembly  two  or  more  times  and  have  the  sampling  pulses  for 
the  matched  filters  staggered  in  time.  For  example,  if  there 
were  three  separate  DMF  assemblies  the  sampling  pulses  to  each 
would  be  offset  from  each  other  by  one  third  of  a  chip  time. 

This  would  ensure  that  at  least  one  DMF  was  being  sampled  with 
a  timing  error  no  greater  than  one-sixth  of  a  chip  period,  i.e. , 
within  one-sixth  of  a  chip  period  of  the  end  of  each  chip. 

3 . 1  Gaussian  Noise 

It  is  shown  in  Appendix  A  (equation  17)  that  when  the  bi-phase 
modulated  signal  is  received  together  with  white  Gaussian  noise  and 
the  signal  is  correctly  aligned  in  the  DMF  the  probability  density 
function  (PDF)  for  the  output  voltage  Z  of  the  DMF  is  given  by 


p(Z) 


Z 

—  exp 
n 


(1) 


where  n  =  the  number  of  chips  in  the  transmitted  sequence 

g 

(■—)  =  signal  to  noise  ratio  in  the  receiver  IF 

N  i 

I  (x)  =  modified  Bessel  function  of  zero  order  and  of  argument  x. 
o 

The  PDF  for  the  output  voltage  of  the  DMF  when  a  sequence  orthogonal 
to  that  to  which  the  DMF  is  matched  is  obtained  by  putting  (S/N)^  =  0 
in  equation  1,  i.e.  , 


p(Z) 


Z 

=  —  exp 
n 


(2) 


This  is  also  the  PDF  for  the  output  voltage  from  the  DMF  for  a  noise 
only  input. 
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3.2  Random  Phase  Sinusoidal  Interference 


It  is  shown  in  Appendix  B  (equation  16)  that  when  the  bi-phase 
modulated  signal  is  received  together  with  an  inband  interfering  sinu¬ 
soid  of  random  phase  and  the  signal  is  correctly  aligned  in  the  DMF, 
the  PDF  for  the  output  voltage  Z  of  the  DMF  is  given  by 


2 

p(Z)  =  —  exp 


1_ 

2n 


(3) 


where  n 


*  number  of  chips  in  the  transmitted  sequence 

=  ratio  of  signal  power  to  interference  power  measured 
in  the  receiver  IF 


I  (x)  =  modified  Bessel  function  of  zero  order  and  argument  x. 

Similarly  for  a  random  phase  sinusoidal  interference  input  only,  or  a 
random  phase  sinusoidal  interference  input  plus  a  bi-phase  modulated 
sequence  orthogonal  to  the  sequence  to  which  the  DMF  is  matched,  the 
PDF  of  the  output  voltage  Z  of  the  DMF  is  given  by 


p(Z) 


z 

—  exp 
n 


(4) 


3 . 3  Fixed  Phase  Sinusoidal  Interference 

It  is  shown  in  Appendix  C  that  for  fixed  phase  sinusoidal  inter¬ 
ference  the  output  of  the  matched  filter  can  be  represented  by  two 
distinct  distributions.  Again  we  assume  that  the  received  sequence 
is  correctly  aligned  in  the  matched  filter.  One  distribution  occurs 
when  the  phase  of  the  interfering  signal  is  such  that  it  causes  both 
the  in-phase  and  quadrature  components  of  the  demodulated  signal  to  be 
completely  suppressed.  We  call  this  effect  double  suppression  and  it 
occurs  when  the  phase  angle  of  the  interfering  sinusoid  x  is  such  that 
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(5) 


X  <  X  <  x± 


where 

-1  A 

X1  =  COS  (—  cos  0)  (6) 

Xq  =  sin  1  (-|  sin  0)  (7) 

and 

A  =  amplitude  of  the  received  signal 
B  =  amplitude  of  the  interfering  sinusoid 
0  =  phase  of  the  received  signal 

When  double  signal  suppression  occurs  the  output  Z  of  the  DMF  is  noise 
only  and  has  a  Rayleigh  PDF  as  given  by  equation  4, 

The  second  distribution  occurs  when  the  phase  of  the  interfering 
signal  is  such  that  it  causes  only  the  inphase  or  the  quadrature  com¬ 
ponent  of  the  received  signal  to  be  completely  suppressed.  This  occurs 
when 


X  <  (In-phase  channel  suppressed) 


X  > 


(Quadrature  channel  suppressed) 


where  x^  and  X^  are  as  previously  defined. 

Equation  19  of  Appendix  C  gives  the  probability  density  function 
for  the  output  Z  of  the  DMF  receiver  when  the  phase  of  the  interfering 
sinusoid  causes  only  one  of  the  quadrature  components  to  be  suppressed, 
i.  e.  , 
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z 


2 


2n 


n 


2 


dZ  (Z  >  n)  (8) 


The  output  Z  can  never  be  less  than  n  since  the  quadrature  component 
which  is  not  suppressed  will  always  give  an  output  of  +  n  which,  when 
combined  with  the  noise  component  from  the  other  quadrature  channel 
(by  taking  the  square  root  of  the  sum  of  the  squares)  will  always  result 
in  Z  >  n. 
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SECTION  IV 


SYNCHRONIZATION  PERFORMANCE 


Each  transmission  in  the  TDMA  system  will  start  with  a  synchron¬ 
ization  preamble  consisting  of  an  n  digit  segment  of  a  pseudo  random 
sequence.  The  n  ones  and  zeros  in  the  synchronization  preamble  bi-phase 
modulate  the  transmitted  carrier  and  are  demodulated  at  the  receiver  by 
a  digital  matched  filter  as  shown  in  Figure  1.  The  particular  n  digits 
of  the  PR  sequence  used  for  synchronization  is  assumed  to  be  stored 
in  the  shift  registers  denoted  by  "code  replica"  in  Figure  1.  Arrival 
of  the  synchronization  preamble  is  assumed  when  the  output  voltage  Z 
of  the  matched  filter  exceeds  a  preset  threshold.  The  system  parameters 
which  determine  just  where  the  threshold  is  set  are  the  probability  of 
detecting  synchronization  P^  and  the  false  alarm  probability  Pp.  These 
two  quantities  are  in  turn  related  to  the  probability  density  function 
(PDF)  of  the  output  voltage  of  the  digital  matched  filter  (DMF).  This 
will  be  shown  in  the  next  three  sections. 


4.1  Gaussian  Noise  Performance 


Equation  1  gives  the  PDF  of  the  output  voltage  Z  of  the  DMF  when 
the  synchronization  sequence  is  correctly  aligned  in  the  DMF.  Equation 
1  may  be  written  as 


p(Z) 


Z 

—  exp 
n 


1_ 

2n 


(9) 


The  probability  of  detecting  synchronization  PD  is  the  probability  that 
Z  exceeds  some  preset  threshold  Z^,  that  is 


P 


D 


oo 


p(Z)  dZ 


(10) 
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Using  equations  9  and  10  we  have 


Making  substitutions , 


I 

o 


dZ  (11) 


(12) 


and 


y-  ’ 


(13) 


equation  11  may  be  rewritten  as 


f 

f  v  exp  - 

i 

+ 

CM 

> 

4 _ 

I 

o 

h  (t,i 

2 

v\ 

-  „ 

dv  (14) 


If  we  let 


o  fS\  ^ 

2n  <K)1  ‘  al 


(15) 


and 


K 


2  =  2 
1  7T 


(16) 


Equation  14  may  be  written  as 


v  exp  - 


r  2  2 

rv  +  (K^)^ 


vKiai 


dv 


(17) 


This  equation  is  the  Marcum  Q  function  and  is  tabulated  in  [5] .  Equa¬ 
tion  17  gives  the  same  probability  of  detection  one  would  obtain  for 
the  envelope  of  a  sinusoid  plus  Gaussian  noise  where  the  input  signal- 


12 


to-noise  ratio  was  given  by 


,S,  2  .Sv 

TT  11  V. 
e  i 


(18) 


With  the  threshold  set  at  6,  the  probability  of  false  alarm  in  the  noise 

2 

only  case  is  given  by  putting  a^  31  0  in  equation  17,  i.e.. 


(19) 


Figure  3  is  a  set  of  receiver  operating  characteristics  (ROC’s)  for  the 
envelope  of  a  sinusoid  in  Gaussian  noise  for  various  values  of  (S/N)^. 
These  ROC!s  are  a  plot  of  the  probability  of  detection  P^  versus  probab¬ 
ility  of  false  alarm  Pp  at  a  fixed  signal  to  noise  ratio  (S/N)e.  For 
convenience  a  threshold  scale  is  provided  along  the  abscissa  and  indicates 
where  the  threshold  at  the  output  of  the  receiver  should  be  set  to  obtain 
specific  values  of  Pj)  and  Pp.  The  units  used  for  the  theshold  are  in 
terms  of  multiples  of  the  rms  noise,  a,  into  the  envelope  detector.  In 
the  next  section  we  illustrate  how  they  can  be  used  to  assess  synchroni¬ 
zation  performance. 

4.1.1  Use  of  the  Receiver  Operating  Characteristics  (ROC)  to 

Assess  Synchronization  Performance  in  Gaussian  Noise 

In  designing  a  synchronization  system  based  on  using  a 
digital  matched  filter  the  quantities  to  be  specified  are: 

a.  the  synchronization  sequence  length,  n 

b.  the  synchronization  input  signal  to  noise 
ratio  (S/N)-£ 

c.  where  the  threshold,  ZT,  should  be  set  at  the 
output  of  the  digital  matched  filter. 
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Figure  3  RECEIVER  OPERATING  CHARACTERISTICS  FOR  THE  DIGITAL  MATCHED 

FILTER  DURING  SYNCHRONIZATION 


These  three  quantities  are  in  turn  dependent  on  the  required  probab¬ 
ilities  of  detection  aid  false  alarm,  P^  and  Pp  which  are  determined 
from  such  operational  requirements  as  message  error  rate,  probability 
of  two  consecutive  messages  being  in  error,  etc. 

We  therefore  assume  that  P^  and  Pp  are  specified  and  that  we 
require  the  three  quantities  given  in  a,  b,  and  c  in  order  to  design 
the  DMF  to  be  used  for  synchronization.  As  a  specific  example  in  use 
of  the  ROC?s  assume  we  require 


P 


D 


.99 


and 


P 


F 


10 


-7 


Looking  at  Figure  3  we  see  that  the  lower  curve  denoted  by 

(S/N)^  =*  15  dB  passes  through  the  intersection  of  the  lines  P^  -  .99 

and  P^  =  10  ^  so  that  we  require  at  least  (S/N)  =  15  dB  to  obtain 
F  e 

the  required  values  of  P^  and  Pp. 


Using  equation  18  we  have 


2 

=  -  n  (tt) 

TT  N 


(20) 


Recalling  that  the  factor  2/tt  represents  a  loss  of  2  dB  we  have  for 

(S/N)  =  15  dB,  that 

e 

n  (|b  =  17  dB  (21) 

N  i 

If  we  now  choose  a  reasonable  length  sequence  n  so  that  the  DMF 
does  not  require  too  many  shift  register  stages,  both  (S/N)^  and  the 
threshold  value  are  fixed. 
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Assume  n  =  128  which  using  equation  21,  gives  the  required  input 
signal  to  noise  ratio 


(f)  -  -* dB 

i 


Looking  at  Figure  3  we  see  that  the  threshold  Z ^  corresponding 


to  P  =  10-7,  is 


ZT  =  5.63  a 


(22) 


Since  a  =  y"n,  the  required  threshold  Z  is  given  by 


ZT  =  5.63  y/l28 
=  63.69 

=  64  (nearest  integer) 


.-7 


In  summary,  if  we  require  P  =  .99  and  P  =  10  and  choose  a  DMF  of 

D  r 

length  n  =  128  we  require  an  input  signal  to  noise  ratio  of  at  least 
-4  dB  and  the  threshold  Z^  should  be  set  at  64. 

4 . 2  Random  Phase  Sinusoidal  Interference 

Equation  3  gives  the  PDF  of  the  output  voltage  Z  of  the  DMF  when 
the  synchronization  sequence  is  correctly  aligned  in  the  DMF  and  the 
bi-phase  signal  is  received  together  with  a  random  phase  sinusoidal 
interfering  signal. 

Equation  3  may  be  written  as 


P  (Z)  =  —  exp  - 
n 


2n 


z2  +  4n 


(S)| 
2  VI 


(23) 


The  probability  of  detecting  synchronization  Pn  is  the  probability 


that  Z  exceeds  some  preset  threshold  Z^,  that  is 
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(24) 


p(Z)  dZ 


Using  equations  23  and  24  we  have 


,.A 
°  Jz  “ 


z 

i 

|*2  ♦*!«§> 

Z 

l~2  <!>  ' 

—  exp  - 
n  r 

2n 

I 

O 

n  \ 

7T 

- 

dZ 

(25) 

We  now  make  the  same  substitutions  that  were  made  in  Section  4.1,  i.e., 
4  =  v  (26) 


=  6 


Equation  25  may  now  be  written  as 


/ 


pD  =7  v  exp  ~ 


rv  +  T 

r 

r  i 

7T 

T 

1T  ^  J 

2 

1 

O 

V 

- 

V 

TT 

Again  if  we  let 


2n  (I)  =  a/ 


and 


K 


2  =  2_ 

2  2 
IT 


equation  28  may  be  written  as 


00 


exp 


V  +  (K2a2)‘ 


I 

vLa. 

O 

2  2 

. 

(27) 


(28) 


(29) 


(30) 


(31) 
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The  function  on  the  right  hand  side  of  equation  31  is  the  Marcum  Q 
function  and  is  tabulated  in  [1].  This  equation  gives  the  same  pro¬ 
bability  of  detection  which  one  would  obtain  for  the  envelope  of  a 
sinusoid  plus  Gaussian  noise  where  the  input  signal-to-noise  ratio  is 
given  by 


(32) 


which  using  equations  29  and  30  gives 


2  n  V 

TT 


(33) 


Again  with  the  threshold  set  at  6  the  probability  of  false  alarm  in 

the  case  where  the  synchronizing  sequence  is  misaligned  in  the  DMF 

and  the  signal  is  received  together  with  the  interfering  signal  is 

2 

obtained  by  putting  a^  =  0  in  equation  31,  i.e.. 


P 


F 


dv 


(34) 


Figure  3  can  be  used  to  assess  synchronization  performance  of  the  DMF 
in  exactly  the  same  way  as  that  used  for  Gaussian  noise  in  Section  4.1.1 
provided  the  different  loss  factor 


K22  =  ^2  =  6.94  dB 

7T 


is  used. 

Using  the  previous  example  all  numbers  would  be  the  same  except 
the  required  input  signal  to  noise  ratio  which  is  obtained  from 


18 


or 


e 


=  15  dB  +  6.94  dB  -  21  dB 


=  +  .94  dB 


In  summary  if  we  require  a  =  .99  and  =  10  ^  and  choose  a  DMF 

D  r 


of  length  128  we  require  an  input  signal  to  interference  power  ratio 
(S/J)  of  at  least  +.94  dB  and  the  threshold  should  be  set  at  64. 

4 . 3  Fixed  Phase  Sinusoidal  Interference 

The  only  case  in  which  synchronization  is  possible  is  when  the 
phase  of  the  interfering  sinusoid  is  in  the  single  signal  suppression 
zone  (i.e.,  only  the  in-phase  or  the  quadrature  component  of  the  signal 
is  suppressed,  not  both).  When  this  occurs  and  the  synchronization 
sequence  is  correctly  aligned  in  the  DMF,  the  PDF  for  the  output  voltage 
Z  is  given  by  equation  8,  i.e., 


(35) 


Z  >  n 


The  probability  of  detecting  synchronization  is  given  by  the  probability 
that  Z  exceeds  some  threshold  voltage  Z T.  That  is 


oo 


dZ 


(36) 


using  equations  35  and  36  we  have 
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p 


D 


dZ 


Making  the  substitution 


Z 


2 


2n 


2 

x 


equation  37  becomes 


(37) 


(38) 


Equation  38  shows  that  if  the  threshold  Z is  set  less  than  n  the 
probability  of  detecting  synchronization  is  unity.  This  point  has  al¬ 
ready  been  covered  in  the  paragraph  immediately  after  equation  8. 


The  PDF  for  the  output  of  the  DMF  when  the  interfering  sinusoid 
or  the  interfering  sinusoid  plus  a  misaligned  synchronization  sequence 
is  received  is  given  by 


p(Z) 


z 

—  exp 


(39) 
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If  we  have  the  same  threshold  voltage  at  the  output  of  the  DMF  the 
probability  of  false  alarm  is  given  by  P^  where 

r 


(40) 


Equation  38  gives  the  probability  of  detecting  the  synchronization 
sequence  given  that  the  phase  angle  of  the  interfering  sinusoid  is 
in  one  of  the  single  signal  suppression  zones.  To  obtain  the  over¬ 
all  probability  of  synchronizing  we  must  weight  the  probability  of 
detecting  the  synchronization  sequence  by  the  probability  that  the 
phase  angle  of  the  interfering  sinusoid  is  in  one  of  the  single  signal 
suppression  zones.  To  do  this  requires  that  we  make  some  assumptions 
regarding  the  distribution  of  the  phase  angle  of  the  interfering  sinu¬ 


soid  . 


A  reasonable  assumption  is  that  the  phase  angle  remains  constant 
for  the  duration  of  a  synchronization  sequence  (less  than  100  usee) 
and  that  the  phase  angle  for  consecutive  synchronization  transmissions 
are  independent  and  are  uniformly  distributed  (0  <  x  <  2tt)  .  Based  on 
this  assumption  the  probability  that  the  phase  angle  lies  in  the  single 
signal  suppression  zone  is  given  by  P^,  where 


P^  =  1  -  (probability  of  double  signal  suppression) 


It  is  shown  in  Appendix  C  (Section  C.3)  that  the  probability  of  double 
signal  suppression  is  given  by 


P  =  —  cos  0) 

S  7T 


(41) 
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where 


cos  w 


(42) 


In  equation  42 


and 


A  =  amplitude  of  the  received  signal 
B  =  amplitude  of  the  interfering  signal 
0  =  phase  of  the  received  signal 
co  =  size  of  the  angular  zone  in  which  the  phase 
of  the  interfering  signal  must  lie  to  cause 
double  suppression 


It  is  also  shown  in  Section  C.3  that  the  angular  zone  co  has  its  maxi- 

mem  size  co  when  0  =  0  or  ti/2  and  its  minimum  size  go  ,  when  0  =  tt/4, 
max  min  5 

therefore  using  equations  42  and  43,  we  have 


COS  CO 


max 


A 

B 


(44) 


and 


cos  co  . 
min 


(45) 
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Since 


_S 

J 


equations  44  and  45  may  be  written  as 


cos  0) 


max 


and 


cos  0)  , 
min 


(46) 


(47) 


Since  the  probability  of  double  suppression  is  a  maximum  when  to  -  u> 

max 

and  a  minimum  when  u)  =  u>  ^  we  may  obtain  upper  and  lower  bounds  on 
the  probability  of  double  suppression  by  using  equations  41  with  46 
and  41  with  47,  that  is 


As  stated  previously,  the  probability  that  the  phase  angle  of  the  inter¬ 
fering  signal  lies  in  the  single  signal  suppression  zone  is  given  by 
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/ 


=  1  -  (probability  of  double  signal  suppression) 


-  1  -  Ps  (50) 

We  may  therefore  obtain  upper  and  lower  bounds  on  P^  by  using  equation 
50  with  48  and  49,  that  is 


P  =  1  -  P 

l,max  s,min 


and 


P 


1  ,min 


1  -  P 

s  ,max 


=  1 


2  -1 
—  cos 

7T 


(51) 


(52) 


Using  the  previous  work  we  may  now  obtain  upper  and  lower  bounds  on 
the  probability  of  synchronizing  since  the  probability  of  synchronizing 
P(sync)  is  given  by  the  probability  of  detecting  sync  given  that  the 
phase  angle  of  the  interfering  sinusoid  is  in  one  of  the  single  signal 
suppression  zones ,  that  is 


P (sync)  =  PD  •  P1 


(53) 


so  that 


P (sync) 


max 


=  P, 


1  ,max 


and 


(54) 


min 


P (sync) 


P  •  P 
D  l,min 


(55) 


Using  equations  54,  51  and  38,  we  have  as  an  upper  bound  on  the  pro¬ 
bability  of  synchronizing 


P  (sync) 


max 


L2  ~  2^ 

^  -  -  V-v1 


1  2  -1 
1 - cos 

TT 


(56) 


Also  using  equations  55,  52  and  38  we  have  as  a  lower  bound  on  the 
probability  of  synchronizing 


P (sync) 


min 


1  -  erf 


/v7"2] 

r 

/  _T_  . 

»  2n 

i  2  -i 

1  -  “  COS  A  /  (7) 
IT  \/  J 


(57) 


The  probability  of  false  alarm  with  the  threshold  set  at  is  still 
given  by  equation  40,  i.e., 

_  2 


pF  =  exp  -  (2^r} 


(58) 


4.3.1  Comments  on  Synchronization  When  the  Interfering  Signal 
is  a  Fixed  Phase  Sinusoid 

Looking  at  the  equations  for  the  probability  of  synchroni¬ 
zation  (equations  56  and  57)  we  see  that  there  is  no  increase  in  the 
probability  of  sync  if  we  set  the  threshold  Z ^  less  than  n  since  the 
first  factor  is  unity  for  Z^  <_  n.  For  Zj  _<  n  the  probability  of  syn¬ 
chronization  is  controlled  solely  by  the  second  factor  which  is  the 
probability  that  the  phase  angle  of  the  interfering  sinusoid  is  in 
one  of  the  single  signal  suppression  zones  and  depends  on  (S/J)  onlv. 

If  we  set  the  threshold  Z^  =  n  we  see  that  using  equation 
58  the  probability  of  false  alarm  is  given  by 


Pp  =  exp  - 


(59) 
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which  equals 


1.2  x  10 


-14 


for  n  =  64 


1.6  x  10 


-28 


for  n  =  128 


and 


2.5  x  10 


-56 


for  n  =  256 


We  see  therefore  that  for  the  threshold  set  at  n  the 
probability  of  false  alarm  is  extremely  small  for  the  values  of  n 
which  would  be  considered  and  that  the  probability  of  synchronization 
is  determined  solely  by  the  second  term  in  equations  56  and  57. 

Figure  4  is  a  plot  of  upper  and  lower  bounds  on  versus  (S/J)  which 
for  _<  n  are  identical  to  the  upper  and  lower  bounds  on  P(sync). 

In  summary,  decreasing  the  threshold  below  Z^  «  n  will  in¬ 
crease  P_,  without  increasing  P(sync).  Increasing  the  threshold  above 
r 

n  will  decrease  P(sync)  and  decrease  P^  which  is  extremely  small  any¬ 
way.  We  therefore  conclude  that  the  threshold  Z^  should  be  set  equal 
to  n  since  this  maximizes  the  probability  of  synchronization  and  results 
in  a  probability  of  false  alarm  P  which  is  extremely  small. 
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Figure  4  UPPER  AND  LOWER  BOUNDS  ON  THE  PROBABILITY  P,  THAT  THE  PHASE 
ANGLE  OF  THE  INTERFERING  SINUSOID  LIES  IN  A  SINGLE  SIGNAL  SUPPRESSION 

ZONE  VERSUS  (y) 
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SECTION  V 


DATA  PERFORMANCE 


In  this  section  we  compute  the  binary  error  rate  for  the  receiver 
described  in  Section  2.2  when  the  received  signal  is  received  together 
with 


a.  white  Gaussian  noise 


b.  random  phase  sinusoidal  interference 

c.  fixed  phase  sinusoidal  interference 


As  described  in  Section  2.2  the  receiver  contains  two  digital  matched 
filters  as  shown  in  Figure  2,  one  filter  being  matched  to  the  MzeroM 


sequence  and  the  other  matched  to  the  "one”  sequence.  If  the  receiver 


is  synchronized  the  outputs  from  the  two  matched  filters  would  be  com¬ 
pared  at  the  end  of  a  received  information  bit  and  a  decision  as  to 
the  reception  of  a  "one"  or  zero  would  be  based  on  the  larger  of  the 
two  DMF  outputs. 

5 . 1  Gaussian  Noise 

When  the  bi-phase  modulated  signal  is  received  together  with 
Gaussian  noise  the  PDF  for  the  filter  output  Z  which  is  matched  to 
the  received  signal  is  given  by  equation  1,  i.e., 


(60) 


The  PDF  for  the  output  voltage  Z-^  of  the  filter  matched  to  a  sequence 
orthogonal  to  the  first  sequence  is  given  by  equation  2,  i.e, 


(61) 
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The  bit  error  probability  is  then  given  by  the  probability  that 
Z^  >  Z  or 


P  =  Prob  (Z,  >  Z) 

e  1 


If  we  let  the  term 


2 

c 


(62) 


(63) 


in  equation  60  we  may  use  the  result  given  by  equation  DIO  of  Appendix  D, 
that  is 

2 

pe  =  7  exp  ”  (t^ 

2 

Substituting  in  equation  64  for  c  from  equation  63  we  obtain 


P 

e 


1 

2 


exp  - 


-  (-> 

TT  V, 


(65) 


So  far  it  has  been  assumed  that  both  the  IF  and  low  pass  filter  band- 
widths  of  the  receiver  are  wide  enough  to  allow  the  pulse  envelope  in 
both  the  I  and  Q  channels  to  come  up  to  full  amplitude  by  the  end  of 
a  chip  time.  If  we  now  wish  to  write  equation  65  in  terms  of  E^/N^  we 
must  make  some  assumptions  regarding  system  bandwidths.  If  we  assume 
that  the  IF  bandwidth  W_^  is  given  by 


W 


i 


(66) 


where  t  is  the  chip  duration  we  have  that  the  energy  per  bit  is 
given  by 


E 


b 


n 


t 
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so  that 


A2  n  .S<. 

N  2N  W, 
o  o  1  i 


Using  equations  65  and  67  we  have 


(67) 


(68) 


For  ideal  quadrature  detection  the  corresponding  expression  for 
probability  of  error  would  be  [2] 


P 


e 


1 

2  exP 


(69) 


Figure  5  shows  a  plot  of  the  probability  of  bit  error  versus  E^/N^ 
for  the  DMF  using  one  bit  quantization.  The  curve  for  ideal  quadra¬ 
ture  detection  is  shown  for  comparison. 

5.2  Random  Phase  Sinusoidal  Interference 


The  PDF  for  the  output  voltage  Z  of  the  filter  matched  to  the 
incoming  signal  in  the  presence  of  random  phase  sinusoidal  interfer¬ 
ence  is  given  by  equation  3,  that  is 


2 

p(Z)  =  —  exp 


1_ 

2n 


(70) 


The  PDF  for  the  output  voltage  of  the  filter  matched  to  a  sequence 
orthogonal  to  the  first  sequence  and  receiving  this  sequence  together 
with  random  phase  sinusoidal  interference  is  given  by  equation  4,  that 
is 

z  z  2 

P(Z1)  =  ^  exp  -  (^-)  (71) 
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Figure  5  BIT  ERROR  PROBABILITY  VERSUS  Eb/N0  FOR  A  DIGITAL 
MATCHED  FILTER  USING  ONE  BIT  QUANTIZATION 
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The  bit  error  probability  is  then  given  by  the  probability  that 
>  Z  or 


P  =  Prob  (Z_  >  Z) 
e  1 

As  in  the  previous  section  if  we  now  let  the  term 


7T 


(72) 


in  equation  70  we  may  again  use  the  result  given  by  equation  DIO  of 
Appendix  D,  that  is 


P 

e 


1  /C  V 

2  exp  "  W 


(73) 


Substituting  in  equation  73  for  c  from  equation  72  we  have 


P 

e 


1 

2  exp 


—  (-) 

2  V 

7T 


(74) 


Figure  6  is  a  plot  of  the  bit  error  rate  P^  versus  the  product  of 
the  processing  gain  n  and  the  signal  to  interference  power  ratio  (S/J). 

5 . 3  Fixed  Phase  Sinusoidal  Interference 

It  has  already  been  shown  in  Appendix  C  and  summarized  in  Section 

3.3  there  are  only  two  possible  states  for  the  output  of  the  digital 
matched  filter  receiver.  The  first  state  occurs  when  the  phase  angle 
of  the  interfering  sinusoid  lies  in  the  double  signal  suppression  zone 
causing  suppression  of  both  the  inphase  and  quadrature  components  of 
the  signal.  For  this  state,  the  output  of  the  receiver  is  noise  only 
resulting  in  a  purely  random  output  bit  stream.  The  probability  of 

a  bit  error  for  this  condition  is  then  given  by 
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Figure  6  BIT  ERROR  PROBABILITY  VERSUS  THE  PRODUCT  OF  PROCESSING 
GAIN  AND  SIGNAL  TO  JAMMING  RATIO  Pe  VERSUS  "  { -j-) 

! 
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(75) 


The  second  state  occurs  when  the  phase  of  the  interfering  sinusoid 
lies  in  one  of  the  single  signal  suppression  zones  causing  suppression 
of  either  the  inphase  or  quadrature  component  of  the  signal  but  not  both. 
For  this  second  state,  the  probability  of  a  bit  error  is  given  by  equa¬ 
tion  (26)  of  Appendix  C,  i.e. , 


(76) 


5.3.1  Comments  on  the  Bit  Error  Rate  Against  Fixed  Phase 
Sinusoidal  Interference 

(a)  For  values  of  n  likely  to  be  used  in  the  system  (n  = 
64,  128,  256)  the  value  of  Pe  as  given  by  equation  (76)  is  extremely 


small  ranging  from  10  to  10  .  It  is  therefore  apparent  that  the 


error  rate  performance  of  the  receiver  would  start  to  be  controlled 
by  the  "front  end"  receiver  noise  well  before  an  error  rate  like  10 
could  ever  be  reached. 


(b)  A  necessary  condition  for  a  usable  bit  error  rate  is 


of  course  that  the  receiver  by  synchronized.  If  the  receiver  does  not 
detect  the  synchronizing  preamble  the  message  will  be  missed  since  the 
timing  circuits  which  control  the  sampling  of  the  data  portion  of  the 
transmission  will  not  be  enabled.  We  saw  in  Section  4.3  particularly 
Figure  13  that  the  probability  of  synchronizing  against  fixed  phase 
sinusoidal  interference  was  low,  having  upper  and  lower  bounds  of 
.52  and  .38  for  a  signal  to  jamming  power  ratio  of  (S/J)  =  -5  dB. 

The  upper  and  lower  bounds  at  (S/J)  =  -10  dB  are  .28  and  .20  respectively. 
We  see  therefore  that  for  this  type  of  interference  at  a  signal  to 
jamming  ratio  of  -10  dB,  on  the  average  we  will  miss  between  7  and  8 
messages  out  of  every  10  transmitted. 
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SECTION  VI 


CONCLUSIONS 


The  analysis  of  the  performance  of  a  digital  matched  filter  (DMF) 
receiver  using  1  bit  quantization  in  both  the  inphase  and  quadrature 
channels  of  the  DMF  when  the  input  signal  is  a  biphase  modulated  signal 
in  the  presence  of : 

1)  Gaussian  noise 

2)  random  phase  sinusoidal  interference 

3)  fixed  phase  sinusoidal  interference 

has  given  the  following  results: 

6 . 1  Gaussian  Noise 

a.  There  is  a  loss  of  1.96  dB  (2/tt)  in  detection  efficiency 
compared  with  ideal  quadrature  detection  (analog  integra¬ 
tion).  This  loss  is  solely  attributable  to  the  use  of  1  bit 
quantization  as  opposed  to  analog  integration. 

b.  The  bit  error  rate  for  a  receiver  using  this  type  of 
digital  matched  filter  is  given  by 

i  V 

7T  N 

o 


c.  The  output  voltage  Z  of  the  DMF  has  a  Rician  probability 
density  function  (PDF)  given  by 

2 


rz2  +  4n 


P (Z)  =  —  exp  - 


(!)  - 

N  i 


2n 


I  - 


when  samples  of  the  received  biphase  signal  are  correctly  aligned 
in  the  filter  (the  "in  sync"  condition). 
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The  output  voltage  Z  of  the  DMF  has  a  Rayleigh  PDF 
given  by 


when  the  samples  are  either  Gaussian  noise  alone,  or 
Gaussian  noise  plus  a  biphase  signal  incorrectly  aligned 
in  the  filter  (the  "out  of  sync"  condition) . 


These  two  facts  allow  the  synchronization  performance  of 
the  DMF  to  be  determined  using  standard  [3,4]  receiver 
operating  characteristics. 


6.2  Random  Phase  Sinusoidal  Interference 

2 

a.  There  is  a  loss  in  detection  efficiency  of  6.94  dB  (2/tt  ) 
compared  with  ideal  quadrature  detection  (analog  integration).  This 
loss  is  solely  attributable  to  the  use  of  one  bit  quantization  as 
opposed  to  analog  integration. 

b.  The  bit  error  rate  for  a  receiver  using  this  type  of  digital 
matched  filter  is  given  by 

2  V 

TT 

where 


n  =  receiver  processing  gain  (i.e.,  number  of  chips/bit) 

S 

(—)  =  signal  to  interference  power  ratio  at  the  input  to 
the  receiver. 

c.  The  output  voltage  Z  of  the  digital  matched  filter  has  a 
Rician  probability  density  function  given  by 
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when  samples  of  the  received  biphase  signal  are  correctly  aligned  in 
the  DMF  (the  "in  sync"  condition). 


The  output  voltage  Z  of  the  DMF  has  a  Rayleigh  PDF  given  by 


when  the  samples  are  random  phase  sinusoid  alone  or  a  random  phase 
sinusoid  plus  a  biphase  signal  incorrectly  aligned  in  the  filter  (the 
"out  of  sync"  condition). 

Again  these  two  facts  allow  the  synchronization  performance  of 
the  DMF  to  be  determined  using  the  same  set  of  receiver  operating 
characteristics  [3,4]. 

6 . 3  Fixed  Phase  Sinusoidal  Interference 

a.  When  the  interference  is  a  fixed  phase  sinusoid  it  is  shown 
that  there  are  only  two  possible  states  for  the  output  of  the  digital 
matched  filter  for  the  case  where  the  received  sequence  is  correctly 
aligned  in  the  matched  filter.  One  state  is  that  in  which  the  phase 
of  the  interfering  sinusoid  causes  both  the  inphase  and  quadrature 
components  of  the  demodulated  signal  to  be  completely  suppressed. 

We  call  this  effect  double  suppression  and  it  occurs  when  the  phase 
angle  of  the  interfering  sinusoid  x  is  such  that 


where 


“I  /A 


(-  cos  e) 


Y .  =  COS 

A1 


X  =  sin  1  cos  6) 


q 
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and 


A  =  amplitude  of  the  received  signal 
B  =  amplitude  of  the  interfering  sinusoid 
0  =  phase  of  the  received  signal 


When  double  signal  suppression  occurs  the  output  of  the  receiver  is 
noise  only.  The  second  state  of  the  receiver  output  is  that  which 
occurs  when  only  the  inphase  or  quadrature  component  of  the  received 
signal  is  suppressed.  This  occurs  when 


X  <  X 
X  >  X 


i 

q 


(inphase  channel  suppressed) 
(quadrature  channel  suppressed) 


where  and  X^  are  as  previously  defined. 

b.  The  bit  error  rate  for  the  digital  matched  filter  receiver 
is  given  by 

Pe  ’I 

when  double  signal  suppression  occurs,  and  by 
pe  yr exp  ■  (f> 

when  single  signal  suppression  occurs. 

c.  When  double  signal  suppression  occurs  it  is  not  possible  for 
the  receiver  to  synchronize  correctly  since  the  receiver  output  is 
noise  only. 

It  is  also  shown  that  the  probability  of  synchronizing  when 
the  samples  of  the  signal  are  correctly  aligned  in  the  DMF  has  upper 
and  lower  bounds  given  by 
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These  results  are  based  on  the  assumptions  that 


1-  the  phase  angle  of  the  interfering  sinusoid  remains 

constant  for  the  duration  of  a  synchronizing  sequence, 

2.  the  phase  angle  of  the  interfering  sinusoid  for  consecu¬ 
tive  synchronization  transmissions  are  independent  and 
are  uniformly  distributed  (0  <  x  <  2u) , 

3.  the  voltage  threshold  at  the  output  of  DMF  is  set  at  a 
value  =  n.  This  value  for  is  optimum  as  discussed 
in  Section  4.3.1. 
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APPENDIX  A 


PROBABILITY  DENSITY  FUNCTION  FOR  THE  DMF  OUTPUT  WHEN  THE 
SIGNAL  IS  RECEIVED  TOGETHER  WITH  WHITE  GAUSSIAN  NOISE 


For  a  bi-phase  modulated  signal  plus  Gaussian  noise  the  probability 

density  functions  of  the  inputs  to  the  I  &  Q  channel  are  both  Gaussian. 

2 

The  I  channel  input  has  a  mean  value  A  cos  0  and  a  variance  a  ,  the 

2  1 

Q  channel  has  a  mean  value  A  sin  0  and  variance  a  .  The  PDF’s  are 

q 

shown  in  Figure  7.  Consider  the  I  channel  only;  the  density  function 
for  the  input  is 

2 


,  N  1  ,x  -  A  cos  0v  „ 

p<lt)  exp  -  <7 - >  d* 


(i) 


The  probability  p^  that  a  digit  is  shifted  into  the  register  with  an 
incorrect  sign  (i.e.,  different  to  that  stored  in  its  matching  position 
in  the  register)  is  simply  the  probability  that  it  has  the  opposite 
sign  to  +A  cos  0  or  to  -A  cos  0,  depending  on  the  phase  of  the  received 
signal) .  Therefore 


This  can  be  shown  to  be  equal  to 


pi  =  2 


,  ,A  cos  0v 

1  -  erf  - ) 


n/3*! 


Assuming  that  the  input  signal  to  noise  ratio  is  small,  i.e., 


<<  1 


2  a. 
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Figure  7a  PROBABILITY  DENSITY  FUNCTION  FOR  THE  INPUT  TO  THE 
I  CHANNEL  A/D  CONVERTER 


OJ 

CO 
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Figure  7b  PROBABILITY  DENSITY  FUNCTION  FOR  THE  INPUT  TO  THE 
Q  CHANNEL  A/D  CONVERTER 
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the  first  term  of  the  series  expansion  for  erf  may  be  used  as  a  good 
approximation ,  i .  e .  , 


2  /A  cos  6X 


(4) 


f—  V  o  /  * 

A  2  °i 


Using  equations  3  and  4  we  have 


Similarly,  it  may  be  shown  that  the  probability  of  an  incorrect  digit 
being  fed  into  the  Q  channel  shift  register  is 


(6) 


Equations  5  and  6  give  the  probabilities  that  a  particular  digit 
(1  or  0)  in  the  I  or  Q  shift  register  do  not  match  that  of  the  stored 
sequence . 

Again  restricting  attention  to  the  I  channel,  the  probability 
that  exactly  k  of  the  n  digits  in  the  shift  register  (for  the  in  sync 
condition)  are  correct  is  given  by  the  binomial  distribution. 


(7) 


The  digital  number  or  voltage  V^.  at  the  output  of  the  I  channel  shift 
register  is  given  by 


Vj  =  Number  of  correct  digits  ~  Number  of 


incorrect  digits 


=  k  -  (n  -  k) 
=  2k  -  n 


I 


(8) 
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The  number  k  is  binoraially  distributed  with  a  mean  value  k  given  by 


k  =  n  (1  -  pi) 


and  variance 

°k2  ■  n  a  -  pi)  pi 

Using  equations  8  and  9  the  mean  value  is  given  by 

V  =  2k  -  n 

=  2n  (1  -  p±)  -  n 
Vx  -  n  (1  -  2P±) 


Substituting  for  p^  from  equation  5  in  equation  11  we  have 


—  2n  A  cos  0  2n  /,S N  . 

v!  7^7  -  j*  V(" \  cos  9 


Also  the  variance  ay  is  given  by 

VI 


2  2 

\  -  (2V 


so  that  using  equations  5  and  10  we  have 


=  4n 


n 


ln  ,  2_  A  cos  9  v  1_  n  2_  A  cos  0, 
2  U  +  /s\72"o~  '  o  U  - 


1  ,,  2  A  cos  6X 

1  w  A 


/  2  A  cos  0N^1 

-  J 


~  n 


for  ,(7 r^}  =  ^  <<  1 


(9) 


(10) 


(11) 


(12) 


(13) 
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Using  a  similar  procedure  it  may  be  shown  that  the  output  from  the 

q  channel  shift  register  has  a  mean  and  variance  given  by 


and 


V 

q 


2n 

7T 


A  sin  9 

T  °q 


n 


(14) 


(15) 


For  n  large  and  p^  approaching  .5  (i.e.,  small  input  signal  to  noise 
ratio)  the  binomial  distribution  for  V^.  (and  V^)  may  be  closely  approx¬ 
imated  by  a  normal  distribution  with  the  same  mean  and  variance. 


The  final  output  Z  from  the  digital  matched  filter  is  obtained  by 
taking  the  square  root  of  the  sum  of  the  squares  of  VT  and  V  . 

1  q 

To  assess  the  performance  of  the  DMF,  we  finally  require  the  pro¬ 
bability  density  function  for  Z. 


It  may  be  shown  [5]  that  the  square  root  of  the  sum  of  the  squares 

of  two  normally  distributed  random  variables  with  identical  variances 

2 

a  and  means  m^  and  m^  is  given  by 


p(Z)  =  ^2 

0 


exp 


1 

_2  .  ,  2  ^  2. 

Z  +  (mi  +  m2  ) 

0  2 

2o 

k 

I 

o 


(16) 


Substituting  in  equation  16  for  the  mean  and  variance  of  the  two  normal 
distributions  from  equations  12,  13,  14  and  15  we  obtain  for  the  PDF 
of  the  output  voltage  Z 


*  The  accuracy  of  this  approximation  is  good  only  in  the  region  of 
the  mean  value  of  the  density  function,  the  accuracy  increasing  as 
p  approached  .5.  Towards  the  tail  of  the  density  function  the 
approximation  is  never  good.  See  for  example,  "Probability  and 
its  Engineering  Uses",  T.  C.  Fry,  Pages  227-234. 
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(17) 


This  is 
DMF  for 


the  required  probability  density  function  for  the  output  of  the 
the  in  sync  condition. 
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APPENDIX  B 


PROBABILITY  DENSITY  FUNCTION  FOR  THE  DMF  OUTPUT  WHEN  THE  SIGNAL 
IS  RECEIVED  TOGETHER  WITH  A  RANDOM  PHASE  INTERFERING  SINUSOID 


Figure  8  shows  a  vector  diagram  of  the  received  signal  of  am¬ 
plitude  A  and  the  interfering  sinusoid  amplitude  B.  The  phase  angle 
of  the  interfering  sinusoid  is  assumed  to  have  a  uniform  distribution 
-7T  <  X  <  -Hr,  that  is  during  the  integration  period  of  n  chips  all 
angles  of  <f>  are  assumed  equally  probable*  Figure  9  shows  the  probab¬ 
ility  density  function  of  the  input  x  to  the  I  channel  A/D  converter 
while  Figure  10  shows  that  of  y  the  input  to  the  Q  channel  A/D  con¬ 
verter  • 


Looking  at  Figure  9  we  see  that  the  probability  p_^  that  the 
output  x  from  the  I  channel  A/D  converter  and  hence  the  probability 
that  the  digit  (one  or  zero)  fed  into  the  shift  register  has  a  sign 
different  to  that  stored  in  the  matching  position  in  the  comparison 
register  is  given  by 


which,  due  to  symmetry  of  the  probability  density  function  simplifies 


or 


1  1  ,  -1 
pi  “  2  “  7  Sin 


(A  cos  6) 
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Figure  8  VECTOR  DIAGRAM  OF  BIPHASE  MODULATED  SIGNAL  A 
AND  INTERFERING  SINUSOID  B 


Figure  9  PROBABILITY  DENSITY  FUNCTION  FOR  THE  INPUT  TO 
THE  I  CHANNEL  A/D  CONVERTER 


Figure  10  PROBABILITY  DENSITY  FUNCTION  FOR  THE  INPUT  TO 
THE  0  CHANNEL  A/D  CONVERTER 
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For  B  >  A  (say  two  or  three  times  as  large) 

=  i  -  i  A  cos  8 
Pi  2  ~  it  B 


(1) 


Similarly  the  probability  of  an  incorrect  digit  being  entered  in  the 
Q  channel  is 


_  JL  _  A  sin  6 
2  tt  B 


(2) 


Restricting  attention  to  the  I  channel,  the  probability  that  exactly 
k  of  the  n  digits  in  the  shift  register  (for  the  in  sync  condition) 
is  given  by  the  binomial  distribution 

p(k)  =  (£)  (1  -  p±)k  Pin_k  (3) 


The  digital  number  of  voltage  at  the  output  of  the  I  channel 
shift  register  is  given  by 

V^.  =  Number  of  correct  digits  —  Number  of 
incorrect  digits 

=  k  -  (n  -  k) 

V  =  2k  -  n  (4) 

The  number  k  is  binomially  distributed  with  a  mean  value  k  given 
by 


k  =  n  (1  -  p±)  (5) 


Approximating  the  sine  of  an  angle  by  the  angle  measured  in  radians 
only  leads  to  an  error  of  approximately  5%  for  an  angle  even  as  large 
as  30°.  This  ensures  that  the  approximation  used  will  always  have  an 
error  less  than  5%  if  B  >  2A. 
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and  variance 


°k  =  n  (1  "  pi)  pi 


(6) 


Using  equations  A  and  5 ,  the  mean  value  V ^  is  given  by 


2k  -  n 

2n  (1  -  p±)  -  n 
n  (1  -  2P±) 


(7) 


Substituting  in  equation  7  for  p^  from  equation  1,  we  have 
—  2n  A  cos  0 


(8) 


Also  the  variance  o  is  given  by 

VI 


V, 


(9) 


Using  equations  6  and  9  we  have 


An  (1  -  p±)  pj 


(10) 


Sub 


stituting  for  from  equation  1  in  equation  10,  we  have 


2  ,  ,1  ,  1  A  cos  0N  ,1  1  A  cos  0N 

°V  '  4n  (2  +  V  - B - >  (2~7 - B - > 


=  n 


_  ,2_  A  cos  0  v ‘ 
1  ^TT  B  ^ 


^  n  for  ~  <<  1 

D 


(ID 
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Similarly,  the  mean  and  variance  of  the  output  voltage  from  the  Q 
channel  shift  register  is  given  by 


—  _  2n  A  sin  0 

V  %  —  - - - 

q  tt  B 


(12) 


and 


aq  *  n 


Following  the  procedure  covered  in  Appendix  A,  and  in  particular 
equations  16  and  17,  we  have  that  the  probability  density  function 
for  the  output  Z  of  the  square  root  circuit  is  given  by 


(13) 


p(Z)  =  -  exp  - 


A.  {z2  +  in!  4 

2"  *2  B2 


i  £  —  h 

o  n  tt  B 


(14) 


Since 


B 


-  =  (— ) 

2  \r 


Signal  Power 


Interference  Power 


(15) 


equation  14  may  be  rewritten  as 


z 

1  {z2  +  Ar2  (s)> 

2n  ^  +  2  V' 

z 

<T> 

=  —  exp  - 
n 

I 

o 

n  > 

IT 

(16) 
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APPENDIX  C 


ANALYSIS  OF  VARIOUS  ASPECTS  OF  DMF  PERFORMANCE 
AGAINST  FIXED  PHASE  SINUSOIDAL  INTERFERENCE 


C. 1  Introduction 

In  this  appendix  we  consider  some  aspects  of  the  performance  of 
the  Digital  Matched  Filter  Receiver  against  fixed  phase  sinusoidal 
interference.  In  some  respects  this  work  is  only  of  theoretical 
interest  since  it  is  unlikely  that  the  receiver  would  be  subjected 
to  this  type  of  interference.  For  the  case  of  an  enemy  jammer  using 
an  inband  jamming  sinusoid  it  would  be  very  difficult  for  the  jammer 
to  remain  exactly  on  frequency  and  in  phase  with  the  communication 
signal  for  a  period  of  time  corresponding  to  the  transmission  of  many 
bits  of  information.  This  inability  to  remain  on  frequency  and  in 
fixed  phase  relation  to  the  communication  signal  is  due  to  frequency 
instabilities  in  the  communication  transmitter  and  receiver  together 
with  a  possible  doppler  frequency  due  to  one  or  both  terminals  being 
airborne.  One  possible  case  of  an  on  frequency  fixed  phase  inter¬ 
fering  signal  is  that  of  a  strong  multipath  reflection  of  the  communi¬ 
cation  signal  into  the  receiving  antenna. 

In  summary,  although  it  appears  that  the  case  of  an  on-frequency 
fixed  phase  interfering  sinusoid  is  unlikely  it  should  be  considered 
since  it  will  represent  an  upper  bound  on  performance  against  a  sinu¬ 
soidal  interfering  signal. 

C2 .  Signal  Suppression  Zones 

In  this  section  we  consider  the  conditions  under  which  it  is 
possible  for  a  fixed  phase  sinusoidal  interfering  signal  to  com¬ 
pletely  suppress  the  inphase  or  quadrature  components  of  the  DMF 
output  or  both. 
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Figure  11  shows  the  vector  diagram  for  the  signal  of  amplitude 
A  and  phase  0  and  the  jamming  sinusoid  amplitude  B  and  phase  y. 

The  input  x  to  the  I  channel  A/D  converter  of  the  DMF  (see  Fig.  l]_) 
will  be 


x  =  +  A  cos  0  +  B  cos  x  (1) 

and  the  input  y  to  the  Q  channel  A/D  converter  will  be 


y  =  +  A  sin  0  +  B  sin  x  (2) 

Considering  equation  1  we  see  that  the  sign  of  output  of  the  I  channel 
A/D  converter  is  completely  determined  by  the  interfering  signal  B  if 

|  B  cos  X  I  |  A  cos  0|  (3) 


and  by  the  communication  signal  A  if 
| A  cos  0 |  >  | B  cos  x I 


(4) 


Let  the  value  of  x  (considering  only  the  first  quadrant  for  the  moment) 
at  which  this  change  over  occurs  be  denoted  by  x^ >  where 


X 


i 


-1 

cos 


(“  cos  0) 


(5) 


Similarly,  there  is  a  critical  value  of  x  at  which  the  sign  of  the  Q 
channel  A/D  converter  output  is  determined  either  by  the  interfering 
signal  or  the  communications  signal.  Using  equation  2  this  angle  is 
given  by 


Xq  =  sin 


-1 


^A  sin  0) 


(6) 
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0  CHANNEL 


Figure  II  VECTOR  DIAGRAM  OF  THE  Bl- PHASE  SIGNAL  A  AND 
THE  INTERFERING  SINUSOID  B 


I  CHANNEL 


Figure  12  SHOWING  HOW  THE  FIRST  QUADRANT  IS  DIVIDED 
INTO  THREE  ANGULAR  ZONES 


Figure  13  SHOWING  HOW  THE  SIGNAL  SUPPRESSION  ZONES  ARE 
REPEATED  IN  EACH  QUADRANT 
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Figure  12  shows  how  and 
angular  zones  such  that  for 


divide  the  first  quadrant  into  three 


X  <  X± 
X  >  xq 


the  I  channel  is  suppressed 
the  Q  channel  is  suppressed 


If  x  falls  in  the  zone 


then  both  the  I  and  Q  channel  components  at  the  output  of  the  DMF  are 
suppressed.  As  the  phase  angle  x  of  the  jamming  sinusoid  varies 


between  0°  and  2 tt  the  suppression  zones  for  the  I  and  Q  channels  are 
repeated  in  each  quadrant  as  shown  in  Figure  13.  We  see  that  depending 


on  the  value  of  x  that  either  the  I  or  Q  channel  output  has  an  output 


n  units  large  and  the  other  a  noise  component 


that 


both  the  I  and  Q  channels  have  a  noise  component  of  rms  value  %/n  with 
no  signal  component.  It  is  also  obvious  that  the  change  from  a  signal 
component  in  the  DMF  output  Z  to  all  noise  is  very  sharp  and  occurs 
when  ever  x  moves  into  the  double  suppression  zone  (i.e.,  I  and  Q  both 
suppressed) . 

We  see  therefore  that  unlike  the  previous  cases  considered  (i.e., 
Gaussian  noise  or  a  random  phase  sinusoid)  in  appendices  A  and  B  it  is 
possible  to  get  complete  signal  suppression  if  the  phase  of  the  inter¬ 
fering  sinusoid  lies  in  the  double  suppression  zone, 

C. 3  Probability  of  Double  Suppression 

We  saw  in  the  previous  section  that  if  the  phase  angle  of  the 
jamming  sinusoid  x  falls  in  the  angular  zone 
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then  both  the  I  and  Q  channel  outputs  are  suppressed  and  the  DMF  out¬ 
put  is  just  noise.  It  is  obvious  that  the  probability  of  the  jamming 
phase  angle  falling  in  or  remaining  in  a  double  signal  suppression 
zone  increases  as  the  angular  zone  between  x  =  and  X  “  increases. 
It  is  therefore  of  interest  to  determine  the  size  of  the  double  suppres¬ 
sion  zone  and  then  relate  it  to  the  probability  of  complete  signal 
suppression. 


Let  to  = 


=  angular  width  of  the  double  suppression  zone. 


From  equations  5  and  6  we  have 


and 


-1 

cos 


(—  cos  0) 
B 


xq  =  sin"1  (f  sin  6) 


(7) 


(8) 


Since 


cos  to  =  cos  (x±  -  X  )  =  cos  x±  cos  x  +  sin  x±  sin  xc 


we  have  using  equations  7  and  8  that 


cos  to  = 


n  2  2~  n  2  t~ 

A  cos  9  VB  -  A  sin  0  VB  -  A  cos  9  A  sin  9 


A 

B 

A 

B 


/  A2  2  / A2  2 

cos  0  i  — 2  s^n  ®)  +  sin  eyu  — 2  cos 


cos  (0  -  a) 


(9) 
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where 


a  = 


If  B  is  considerably  larger  than  A,  as  would  be  the  case  for  effective 

77 

jamming  a  ■+  -r  so  that  using  equation  9  we  have 


COS  0) 


~  cos  (e  -  j) 


(10) 


~ k  -  II 

0  <  0  <  2 


Since 


A  _  Signal  Power  _  , Sv 

2  Jamming  Power 

D 

Equation  10  may  be  written  as 


COS  0) 


cos  (0  -  *1*) 


(11) 


For  fixed  (y)  the  double  suppression  zone  has  its  minimum  size 

TT  TT 

when  the  phase  angle  0  of  the  received  signal  is  (y)  ,  i.e.,  cos  (0  -  y)  ■ 
1  so  that  cos  cj  has  its  maximum  value  and  hence  u  its  minimum  value*  The 
suppression  zone  has  its  maximum  size  when  0  5  0  or  y 


If  we  assume  that  the  phase  of  the  jamming  sinusoid  has  a  uniform 
distribution  bit  to  bit  and  remains  essentially  constant  during  an  in¬ 
formation  bit  the  probability  of  complete  signal  suppression  is  given 
by  the  size  of  the  angular  suppression  zone  divided  by  i.e., 


Prob  (Signal  Suppression)  =  Ps 


0) 

tt/2 


57 


using  equation  11  we  have 


cos 


-1 


2  Oj)  cos  (0  -  j) 


(12) 


The  cosine  term  in  equation  12  always  has  a  value  between  .7071  and  1.0  so 
that  a  reasonable  approximation  for  Ps  is 


(13) 


for  (j)  «  1 


C. 4  Probability  Density  Function  for  the  DMF  Output  Voltage  (Single 

Suppression  Zone) 

In  this  section  we  determine  the  PDF  for  the  output  voltage  from 
the  DMF  when  the  phase  angle  of  the  jamming  sinusoid  is  in  one  of  the 
single  suppression  zones. 

The  probability  density  function  for  the  output  Z  of  the  square  root 
circuit  of  the  receiver  (see  Fig.  2)  when  x  is  not  in  the  double  sup¬ 
pression  zone  is  given  by  the  envelope  of  a  constant  amplitude  inphase 
—  2 

component  (V  =  n  and  ay  =  0)  plus  a  quadrature  Gaussian  noise  compon- 

—  ^2  i 

ent  (V  =0  and  av  =  n) . 
q  vq 

Figure  14  shows  a  sketch  of  the  normal  PDF  for  the  quadrature 
output  V  from  the  DMF  located  in  the  plane  V  =  n.  The  probability 

q  i 

that  the  output  Z  of  the  DMF  has  a  value  between  Z  and  Z  +  dZ  is  equal 
to  the  shaded  area  under  the  normal  distribution  as  shown.  That  is 


P (Z)  dZ  =  2p  (Vq)  dVq 


or 


dV 


p(Z)  =  2p 


q 


(14) 
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IB-37,721 


Figure  14  SKETCH  OF  NORMAL  PROBABILITY  DENSITY  FUNCTION  RESULTING 
FROM  ANGLE  x  IN  ONE  OF  THE  SINGLE  SUPPRESSION  ZONES 


Figure  15  PROBABILITY  DENSITY  FUNCTION  FOR  OUTPUT  VOLTAGE  Z 
FROM  DM F  FOR  JAMMER  PHASE  IN  ONE  OF  THE  SINGLE  SUPPRESSION  ZONES 
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2  2  2 

Since  Z  =  V  +  V,  ,  and  =  n,  we  have 
q  i  i 


2  2  2 
Z  «  V  +  n 

q 


Differentiating  equation  15  we  have 


(15) 


ZdZ  =  V  dV 

q  q 


(16) 


or 


Z_ 

V 

q 


Substituting  for  V  from  equation  15  we  have 


(17) 


and 


=  n 


we  have 


p(V 


(18) 
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Using  equations  14,  15,  17  and  18  we  finally  obtain 


p(Z) 


exp 


or 


p(Z)dZ  = 


2 

y2™ 


exp 


(19) 


Equation  19  gives  the  probability  density  function  for  the  DMF  output 
voltage  when  the  phase  angle  of  the  jamming  sinusoid  is  in  one  of  the 
single  signal  suppression  zones  and  is  sketched  in  Figure  15. 

C. 5  Cumulative  Distribution 

In  work  on  synchronization  it  will  be  necessary  to  have  the  cumu¬ 
lative  distribution  for  the  DMF  output.  Using  equation  19  we  have 


Let 


Z 


2 


2n 


2 

x 


then  ZdZ  =  2nx  dx  so  that  equation  20  becomes 
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P<Z„) 


2 


-  X 


dx 


(21) 


C. 6  Bit  Error  Probability  (Single  Suppression  Zone) 

For  binary  transmission  two  orthogonal  PN  sequences  are  used. 

The  receiver  would  contain  two  digital  matched  filters,  one  matched  to 
each  of  the  two  possible  sequences.  A  bit  error  would  occur  if  the 
voltage  from  the  DMF  to  which  the  received  sequence  was  not  matched 
exceeded  that  from  the  matched  filter. 

The  PDF  for  the  output  voltage  of  the  DMF  matched  to  the  incoming 
sequence  is  given  by  equation  19,  i.e., 


p  (Z)  dZ  —  i—rj 


exp  - 


'2  21 
Z  —  n 


2n 


7TT 7 


dZ 


(22) 


The  PDF  for  the  output  voltage  Z^  of  the  DMF  orthogonal  to  the  re¬ 
ceived  sequence  is  given  by  the  Rayleigh  distribution,  i.e., 


P 


<v  ■  ;r exi> 


(23) 
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The  probability  that  exceeds  the  output  Z  from  the  DMF 
to  the  received  sequence  is  given  by 


p^  =  Prob  (Z^  >  Z) 


OO  00 


The  inner  integral  is  given  by 


so  that  equation  24  reduces  to 


00 

,/ 

J  Z=n 


P  (Z) 


exp 


-Z 


2-1 


2n 


dZ 


Substituting  for  p(Z)  from  equation  22  we  obtain 


Let 


so  that 

ZdZ  =  2nx  dx 


matched 


(24) 


dZ  (25) 
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Equation  25  may  now  be  written  as 


00 


From  equation  26  we  see  that  for  reasonably  large  processing 
gains  (i.e.,  n  >  50)  the  error  rate  due  to  the  phase  angle  of  the 
jamming  sinusoid  being  in  one  of  the  single  suppression  zones  is 
small.  The  controlling  term  in  the  error  rate  for  this  case  would 
be  receiver  noise  and  not  the  jamming  signal. 


(26) 
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APPENDIX  D 


EVALUATION  OF  AN  INTEGRAL  RELATED  TO  ERROR  RATE  CALCULATIONS 


In  determining  the  binary  error  rate  for  a  matched  filter 
receiving  a  bi-phase  signal  together  with  Gaussian  noise  or  random 
phase  sinusoidal  interference,  evaluation  of  the  probability  that  one 
random  variable  exceeds  a  second  random  variable  Z  is  necessary. 
The  two  random  variables  have  probability  density  functions  of  the 
form 


p(Z1)  =  —  exp 


(Dl) 


and 


p(Z) 


z 

—  exp 
n 


2  2 
Z  +  c 


2n 


,Z  \ 
(n  C) 


(D2) 


where  c  is  a  known  constant. 
We  therefore  require 


P  =  Prob  (Z-  >  Z) 
e  i 


We  have  therefore  that 


dZ 


(D3) 


Evaluating  the  integral  containing  only,  we  have,  using  equation  Dl  that 


(D4) 
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Substituting  in  equation  D3  for  the  value  of  the  inner  integral  we 
obtain 


P 

e 


[  2Z2  +  c 
[  2" 


l 

o 


,z  \ 

(n  C) 


dZ 


1 

n 


2 


c 


n 


I 

o 


(f  c)  dZ 


(D5) 


(D6) 


An  integral  similar  to  that  in  equation  D6  is  evaluated  in  [6]  that 
is 


/ 


Z  exp  (-p2Z2)  J  (aZ)  — exp  (— %r) 
°  2pZ  4pZ 


(D7) 


Since 


J  [(/^L  a)  Z]  =  I  (aZ) 
o  v  o 


(D8) 


we  have 


/ 


Z  exp  (-p2Z2)  I  (aZ)  dZ  =  exp  (^-y) 

°  4p 


(D9) 


Using  the  result  in  equation  D9  and  letting 


2  1  ,  c 

p  =  —  and  a  =  — 

K  n  n 


equation  D6  reduces  to 


Pe  =  Y  eXP  (_ 


(DIO) 


66 


REFERENCES 


1.  J.  I.  Marcum,  "Table  of  Q-Functions,  "  Rand  Corporation  Research 
Memorandum  339,  AD116551. 

2.  G.  N.  Watson,  Theory  of  Bessel  Functions,  Cambridge  University 
Press  (1944),  page  394. 

3.  G.  H.  Robertson,  "Operating  Characteristics  for  a  Linear  Detector 
of  CW  Signals  in  a  Narrow-Band  Gaussian  Noise,  "  Bell  System 
Technical  Journal,  Vol.  46,  No.  4  (April  1967),  pages  755-774. 

4.  M.  I.  Skolnik,  Introduction  to  Radar  Systems,  New  York,  McGraw- 
Hill  (1962),  page  34. 

5.  A.  Papoulis,  Probability,  Random  Variables  and  Stochastic  Processes, 
New  York,  McGraw-Hill  (1965),  page  196. 


67 


Security  Classification 


DOCUMENT  CONTROL  DATA  -  R  &  D 


fSfruf/fy  classification  of  title,  body  ol  abstract  and  indexing  annotation  must  be  entered  when  the  overall  report  la  c/as 


1  ORIGINATING  ACTIVITY  (Corporate  author) 

2..  REPORT  SECURITY  C  L  A  SSI  P  1  C  *  Tl  ON 

The  MITRE  Corporation 

UNCLASSIFIED 

P.  O.  Box  208 

2b.  GROUP 

Bedford,  Mass.  01730 

3  REPORT  TITLE 


INTERFERENCE  EFFECTS  IN  DIGITAL  MATCHED  FILTERS 

4  OESCRIPTIVE  NOTES  (Type  ol  report  end  inclusive  dates) 

5  authoriSi  ('First  name,  middle  initial,  last  name) 

P.  P.  Bratt 


ft  REPOR  T  D  A  TE 

7a.  TOTAL  NO.  OF  PACES 

lb.  NO  OF  REFS 

DECEMBER  1972 

75 

6 

6a  CONTRACT  OR  GRANT  NO. 

90.  ORIGINATOR'S  REPORT  NUMBER(S) 

F19628-71-C-0002 

b.  PROJECT  NO 

ESD-TR-72-322 

511A 

C. 

9b.  OTHER  REPORT  NO<5)  (Any  other  numbers  that  may  be  assigned 
this  report) 

d 

MTR-2469 

10.  DISTRIBUTION  STATEMENT 


Approved  for  public  release;  distribution  unlimited. 


11  SUPPLEMENTARY  NOTES 

12  SPONSO  RING  Ml  LI  T  AR  Y  ACTIVITY 

Deputy  for  Planning  and  Technology- 

Electronic  Systems  Division  (AFSC) 

L.  G.  Hanscom  Field,  Bedford,  Mass. 

13  ABSTRACT 

This  MTR  contains  a  summary  of  some  of  the  more  important  theoretical  results  related  to  the 
performance  of  a  digital  matched  filter  receiver  designed  to  receive  a  biphase  modulated  signal 
in  the  presence  of  interference.  Both  the  synchronization  and  data  performance  of  the  receiver 
are  evaluated  when  the  signal  is  received  in  the  presence  of  a)  Gaussian  noise,  b)  Random 
phase  sinusoidal  interference,  c)  Fixed  phase  sinusoidal  interference.  Specific  attention  is  pah 
to  the  effect  of  processing  the  received  signal  using  one  bit  quantization. 


DD  ,"“,1473 


Security  Classification 


Security  Classification 


1  4. 

KEY  WO  ROS 

LINK  A 

LINK  B 

LINK  C 

ROLE 

WT 

ROLE 

W  T 

ROLE 

WT 

DIGITAL  FILTERS 

INTERFERENCE  FILTERS 

SIGNAL  PROCESSING 

Security  Classification 


